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Abstract 


Hypersonic boundary layer flows over a circular cone at moderate incidence can support strong crossflow 
instability. Due to more efficient excitation of stationary crossflow vortices by surface roughness, such 
boundary layer flows may transition to turbulence via rapid amplification of the high-frequency secondary 
instabilities of finite amplitude stationary crossflow vortices. The amplification characteristics of these 
secondary instabilities are investigated for crossflow vortices generated by an azimuthally periodic array of 
roughness elements over a 7-degree half-angle circular cone in a Mach 6 free stream. Depending on the local 
amplitude of the stationary crossflow mode, the most unstable secondary disturbances either originate from 
the second (i.e.. Mack) mode instabilities of the unperturbed boundary layer or correspond to genuine 
secondary instabilities that reduce to stable disturbances at sufficiently small amplitudes of the stationary 
crossflow vortex. The predicted frequencies of dominant secondary disturbances are similar to those 
measured during wind tunnel experiments at Purdue University and the Technical University of 
Braunschweig, Germany. 


1 



Introduction 


By virtue of its influence on surface heat transfer, skin friction, and separation characteristics, boundary layer 
transition from laminar to turbulent state can have a major impact on the design and performance of 
hypersonic flight vehicles. Transition over realistic vehicle surfaces is often caused by surface roughness of 
sufficiently large amplitude. However, when the surface is relatively smooth, the laminar-turbulent transition 
process is initiated by linear instabilities of the laminar boundary layer. Typically, the second (or Mack) 
mode instability dominates transition in 2D/axisymmetric boundary layers at hypersonic speeds [1], although 
centrifugal (i.e., Gortler) instabilities may also come into play when the body surface has concave curvature 
along the streamwise direction [2]. Three-dimensional boundary layers involve the additional mechanisms of 
stationary and traveling modes of crossflow instability [3] and, depending on the geometric configuration, the 
attachment line instability as well. 

Regardless of the speed regime, linear stability correlations have been quite successful in predicting the onset 
of transition when a single instability mechanism dominates the transition process. Mixed mode transition 
has been more difficult to predict. Furthermore, measurements of crossflow instability in low-speed boundary 
layers have exposed the shortcomings of applying purely linear predictive models to transition in 3D 
boundary layers by revealing the importance of nonlinear effects during crossflow dominated transition [3]. 
The canonical configuration of a circular cone at angle of attack (AO A) includes the necessary elements to 
study both mixed mode transition and crossflow development in the context of both supersonic [4-12] and 
hypersonic [13-29] boundary-layer flows. Instability mechanisms for the boundary layer flow over an elliptic 
cone have been investigated in Refs. [19] and [29]-[46]. While the elliptic cone configuration supports 
similar instability mechanisms as the circular cone at an angle of incidence, the details of transition patterns 
can be different as noted in Ref [41]. 

The availability of low-disturbance (i.e., quiet) hypersonic facilities opens up the possibility of a detailed 
study of the above mentioned instability phenomena in a controlled environment [21-25]. However, because 
of the complexity of the crossflow transition process, the interpretation of wind tunnel measurements must be 
guided by a detailed computational analysis of the flow field and its stability characteristics. The objective of 
this ongoing effort is to help provide the required characterization in support of wind tunnel experiments so as 
to eventually resolve some of the outstanding issues related to transition mechanisms in three-dimensional 
hypersonic boundary layers over nominally smooth surfaces. 

The development of stationary and traveling crossflow instabilities on circular and elliptic cone configurations 
at a freestream Mach number of 6.0 has been addressed in Refs. [ 1 7-29] . Limited measurements of traveling 
crossflow instability by Borg et al. [47] have yielded encouraging comparison with the linear stability 
predictions by Li et al. [19]. More recent quiet tunnel experiments indicate that similar comparisons may 
soon be possible for the stationary crossflow modes as well. In particular, experimental measurements of 
crossflow instability over a yawed circular cone in two different Mach 6 quiet wind tunnels have been 
reported recently. Each of these investigations is rather remarkable, albeit for different reasons. Ward et al. 
[21] performed measurements using temperature sensitive paint (to measure time averaged surface heat 
transfer) and dynamic surface pressure transducers at a variety of flow conditions (a =3 deg and 4 deg. Re = 
9.84x 10^ to 12.01x10^, both with and without surface roughness, and a range of axial and azimuthal locations 
of unsteady pressure sensors) to provide strong evidence of secondary instability. Craig and Saric [25] 
performed in-depth hot wire measurements of the boundary-layer flow over a selected range of azimuthal 
locations at a fixed freestream condition. The unsteady pressure spectra measured by Ward et al. [21] in the 
BAMQT facility at Purdue University showed two major peaks at most of the measurement locations and 
flow conditions with a = 3 deg and a = 4 deg. The lower- frequency peak near 40 kHz to 50 kHz was 
attributed to traveling crossflow modes and the higher frequency peak ranging from 125 kHz to 500 kHz (i.e., 
nearly three to ten times the frequency of the lower peak) was attributed to the secondary instability of 
stationary crossflow modes. A limited set of wind tunnel runs showed an additional, intermediate-frequency 
peak near 150 kHz, but the character and the origin of the fluctuations associated with this peak could not be 
established. Subsequent measurements by McKieman et al. [22] and Edelman and Schneider [23] extended 
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the above measurements to a = 6 deg. For a narrower set of sensor locations, the peak frequency of the 
measured surface pressure spectra was between 300 and 335 kHz. The peak frequency decreased slightly as 
the sensor location shifted downstream. 

The hot wire measurements by Craig and Saric [24, 25] in the NASA Langley Mach 6 Quiet Tunnel at Texas 
A & M University provided quantitative maps of steady and unsteady mass-flux distributions across the cross- 
section of the stationary crossflow vortices at selected axial locations within the range ofX= 0.35 and 0.4 
meters and azimuthal locations of (j) between 112 and 124 degrees relative to the windward plane of 
symmetry. Similar to the Purdue measurements [21-23], they observed two peaks in the frequency spectra. 
The low frequency peak over 15 to 60 kHz was approximately consistent with the surface pressure spectra of 
Ward et al. [21]. However, the high-frequency peak in the unsteady hot wire measurements involved 
significantly lower frequencies of 80 to 130 kHz. The high-frequency disturbances were concentrated over 
the inclined shear layer that is visible from the overturning stationary mass-flux contours produced by the 
finite amplitude stationary crossflow vortices. The measured mode shapes were remarkably similar to the 
secondary instability mode driven by the azimuthal shear induced by the stationary vortex. According to 
Craig and Saric [25], the frequency of the azimuthal-shear mode of secondary instability scales like 
2Uel5 , which is equal to 120 kHz for the conditions of their 

measurement. Thus, the frequency range of measured high-frequency disturbances in their experiment was in 
line with this estimate. Computations by Oliviero et al. [26] indicated qualitative similarities with the 
stationary crossflow modes measured in Ref [25]; however, significant quantitative differences were apparent 
in the disturbance amplitudes as well as mass flux contours. 

Conventional facility measurements analogous to the quiet-tunnel measurements at Purdue University have 
been reported by Munoz et al. [20]. They used infrared measurements of surface temperature to infer heat 
transfer and surface pressure measurements to infer the frequencies, amplitudes, and wave angles of 
instability wave packets over a Mach 6 yawed circular cone tested at the Technical University of 
Braunschweig. The average spatial extensions of these packets were also mapped out using circumferential 
cross-correlations. For sensors near (j)=90 deg from the windward ray and ^locations between 0.21 m and 
0.36 m, the low frequency lobe was spread over the range between 20 and 50 kHz and the high-frequency 
peak was centered at approximately 300 kHz to 350 kHz. These frequency values are remarkably similar to 
those measured under quiet conditions and, therefore, even though Munoz et al. attributed the high-frequency 
peak to possible second mode instability, it is possible that this peak was actually associated with secondary 
instabilities of the distorted mean flow in the presence of strong stationary crossflow vortices. At the larger 
azimuthal angle of 1 10 deg, Munoz et al. noted higher disturbance amplitudes and slightly reduced center 
frequency of the high-frequency lobe, analogous to the measurements of Ward et al. [21]. Nonlinear 
interactions between the low-frequency disturbances (traveling crossflow modes) and the high-frequency 
disturbances (attributed to second mode perturbations) were also inferred at the most downstream 
measurement location {X= 0.36 m) by using bicoherence analysis. Analysis of cross-correlations between 
multiple pressure transducers indicated that the dominant propagation angle associated with the low- 
frequency disturbances was between 60 and 70 degrees relative to the flow velocity at the edge of the 
boundary layer. 

Computations of secondary instability of crossflow modes in a hypersonic boundary layer were first reported 
by Li et al. [29]. The predicted secondary instabilities with the highest growth rates had frequencies that were 
comparable to those measured in the Purdue experiments [21-23], but the computations also revealed the 
existence of less unstable modes with lower disturbance frequencies that were comparable to those measured 
by Craig and Saric [25]. The unstable secondary modes were shown to have an intricate variety of modal 
structure, with fluctuations concentrated in various subregions of high shear associated with the stationary 
crossflow structures. However, in general, no dominant modes with spatial concentration on the top of the 
inclined shear layer (the mode shape typically associated with the wall-normal shear of the basic state) were 
found. However, the analysis of Li et al. [29] was limited to local secondary instability analyses at selected 
stations and, hence, did not provide any information about the integrated amplification of the secondary 
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disturbances. The present paper rectifies that shortcoming by tracking the evolution of unstable secondary 
modes over the region of interest. 

The remaining sections of this document are laid out as follows. A brief summary of the flow configuration 
of interest is given in Section 1, which also describes the analysis codes used in this study. Nonlinear 
evolution of stationary crossflow modes using direct numerical simulations is described in Section 2. 1 and the 
secondary instability of these finite amplitude crossflow vortices is investigated in Section 2.2. Summary and 
concluding remarks are presented in Section 3. 
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Nomenclature 

/ = frequency of instability waves [kHz] 

Me = boundary-layer edge Mach number 

Moo = freestream Mach number 

n = azimuthal wavenumber, i.e., number of waves across circumference 
U = streamwise velocity [m/s] 

Ue = boundary-layer edge velocity [m/s] 

Uoo = freestream velocity [m/s] 

P OO — freestream pressure [Pa] 

r.m.s. = root mean square 

Re = freestream unit Reynolds number [m'^] 

5 = surface distance along cone generatrix [m] 

t = time [sec] 

T^, adb ~ adiabatic wall temperature [K] 

= wall temperature [K] 

Too = freestream temperature [K] 

N = N-factor of linear instabilities 

X = axial coordinate [m] 

Y = wall-normal distance [m] 

a = angle of attack [deg] 

= azimuthal wavelength [deg or m, depending on context] 

(j) = azimuthal coordinate with respect to windward meridian [deg] 

T = total wall shear stress [Pa] 

TO = azimuthal wall shear stress [Pa] 
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1. Flow Configuration and Analysis Codes 

The flow configuration of interest corresponds to a circular cone at angle of attack to a uniform, hypersonic 
free stream. The primary configuration for the Purdue experiments corresponds to a 0.457 meter long (18 
inches), nominally sharp nosed cone (modeled with a nose radius of lO "^ times the cone length) with a half 
angle of 7 degrees and an angle of attack equal to 6 degrees. The freestream conditions correspond to a free- 
stream Mach number ofM= 6, unit Reynolds numbering = 10.8x 10^ per meter, and free-stream temperature 

Too = 52.44 K. The temperature of the model surface, is equal to 300 K. 

The unperturbed boundary layer flow over the cone is computed using the VULCAN (Viscous Upwind 
aLgorithm for Complex flow ANalysis) software [48]. The code solves the unsteady, conservation equations 
appropriate for laminar or turbulent flow of calorically or thermally perfect gases with a spatially second 
order accurate cell-centered finite volume scheme. A variety of upwind schemes are available for fine-tuning 
of numerical dissipation and shock capturing in order to ensure solution accuracy for stability computations. 
In the present computations, the inviscid fluxes were constructed using the MUSCL scheme, the van 
Albada gradient limiter [49] and the Low Dissipation Flux Split Scheme (LDFSS) of Edwards [50]. The cell 
face gradients required to construct the viscous fluxes were obtained using an auxiliary control volume 
approach that results in a compact viscous stencil that produces a second order accurate approximation of the 
full Navier-Stokes viscous fluxes. The solutions were relaxed in pseudo time to steady state using the 3-D 
ILU(O) scheme [51] with a constant CFL number on the order of 50. Grids of 345x257x353 and 
345x129x353 in the axial, azimuthal, and wall-normal directions, respectively, were used. Exploiting the 
azimuthal symmetry of the boundary layer about the windward and leeward planes, only one half of the 
circular cone was modeled in the boundary layer computation. In previous work, similar computations of the 
mean flow over an elliptic cone configuration were cross-validated against the solutions obtained with the 
LAURA code [52] for various grid sizes [20]. Thus, no extensive grid convergence analysis was deemed 
necessary. As an illustration, however. Fig. 1 shows the boundary layer profiles obtained using two different 
grids. An excellent agreement may be observed between the two sets of profiles within the region of interest. 
Further details on the mean flow solution and its linear stability characteristics may be found in Ref [17]. 



Fig. 1. Effect of grid resolution on mean velocity profiles along (/> = 120 deg. 

The evolution of stationary crossflow disturbances in the boundary layer flow is computed by solving the full 
three-dimensional compressible Navier-Stokes equations in conservation form. To provide a controlled 
setting for the study of secondary instability, the stationary crossflow vortices are excited via azimuthally 
periodic forcing with a single azimuthal wavenumber of ^ = 60. The forcing is introduced via an array of 
localized roughness elements centered atX= 0.18 m. The streamwise shape of the roughness elements 
corresponds to a half-wavelength cosine shape with an axial wavelength of 0.008 m and a peak height 
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perturbation of 5 |Lim. The resulting disturbance field is similar to that observed in the experiments by Craig 
and Saric [25]; however, because the disturbance field in their experiments was excited by uncontrolled 
background disturbances, no attempt has been made to achieve a precise match with their data. The primary 
focus of the present study is to characterize the general behavior of secondary instabilities of stationary 
crossflow vortices in a hypersonic boundary layer, and for that purpose, it is sufficient that the crossflow 
vortices approximate the strength of nonlinear crossflow vortices as gauged by the rollup of mass-flux 
contours across the cross-section and the spatial region (axial and azimuthal locations) of vortices with 
sufficiently large amplitudes. 

The working fluid is assumed to be perfect gas (air) and the usual constitutive relations for a Newtonian fluid 
are used: the viscous stress tensor is linearly related to the rate-of-strain tensor, and the heat flux vector is 
linearly related to the temperature gradient through the Fourier's law. The coefficient of viscosity is computed 
from the Sutherlands's law, and the coeffcient of thermal conductivity is computed by assuming a constant 
Prandtl number Pr = 0.7 1 . A detailed description of the governing equations and their numerical solution is 
given by Wu et al. [53]. The inviscid fluxes from the governing equations are computed using a seventh-order 
weighted essentially non-oscillatory finite-difference WENO introduced by Jiang and Shu [54], the present 
scheme is optimized by means of limiters [53, 55] to reduce the numerical dissipation. Both an absolute 
limiter on the WENO smoothness measurement and a relative limiter on the total variation are employed 
simultaneously during the simulation. The viscous fluxes are discretized using a fourth-order central 
difference scheme and time integration is performed using a third-order low-storage Runge-Kutta scheme 
[56]. The numerical code has been previously applied to simulation of turbulence in hypersonic boundary 
layers [57] as well as laminar- turbulent transition due to crossflow instability in swept wing boundary layers 
[58, 59]. The lessons learned from these simulations were applied to develop the computational grid for 
stationary crossflow evolution in the present work. 

The primary instability computation is performed on a grid consisting of 600x1200x160 points in the axial, 
azimuthal, and radial directions, respectively. Without any loss of generality, the stationary crossflow 
disturbances are assumed to be symmetric about the windward and leeward planes, so that only one half of the 
circular cone was included in the computational domain. The grid points are spaced uniformly in the 
azimuthal direction. 

The secondary instability of the finite amplitude stationary crossflow vortices computed using DNS is 
investigated using two different methods. The first method corresponds to a quasi-parallel analysis of spatial 
instability as described by Li and Choudhari [60, 61], whereas the second technique accounts for the leading- 
order non-parallel effects using the framework of parabolized stability equations as described by Paredes et al. 
[62, 63]. The results presented in this paper are exclusively based on the quasi-parallel method. Predictions 
based on the PSE are deferred to a follow-on paper. 

In principle, the quasi-parallel stability formulation can be extended to include the effects of surface 
curvature. In the present paper, however, those effects are assumed to be a higher order correction that is 
comparable to the effects of mean-flow non-parallelism. Because the wavelength of the crossflow instability 
is comparable to the thickness of the boundary layer, the above assumption is consistent with the high 
Reynolds number asymptotic limit underlying the boundary layer approximation. Instability analysis in the 
absence of surface curvature amounts to unrolling the cone surface onto a plane. If one were to coat the cone 
with paint and roll it on a fiat surface, the resulting printed area in the shape of a fan would represent the 
cone-to-plane mapping (Fig. 2(a)), with a one-to-one correspondence between the points on the cone and 
those on the unrolled fiat surface. The left half of Fig. 2(b) shows the axial velocity contours in a subset of a 
typical cross-section of the cone in the presence of stationary crossflow vortices, and its map onto a domain 
resting on the fiat, unrolled surface is shown in the right half of Fig. 2(b). For later reference, the vortices 
most easily discernible from their surface signature are numbered 1 through 12 in the figure. The azimuthal 
extent of each vortex at a given axial station is determined by the distance between the two points on either 
side of the vortex where the wall shear in the meridian direction is a local maximum. The vortex is then made 
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spanwise periodic via detrending and Fourier low-pass filtering to facilitate local secondary instability 
computations. In the computational volume with the flattened (i.e., unrolled) base, theX-direction is taken as 
the meridian direction, i.e., the direction along any line joining a point on the flat surface and the virtual 
vertex of the fan. The local spanwise direction at a point is along the azimuthal coordinate at a fixed axial 
location over the surface of the cone. The mean flow resulting from these modifications is now amenable to 
both temporal and spatial stability analyses using the procedure described in Refs. [60, 61] in the context of 
crossflow modes in swept- wing boundary layers. 

Because the mean flow modified by the presence of crossflow vortices involves comparable length scales in 
both azimuthal and radial directions (see Fig. 2(b)), the basic state for secondary instability is inhomogeneous 
in two spatial directions. Hence, predicting the amplification characteristics of secondary instability modes 
leads to a planar (i.e., 2D), partial differential equation based eigenvalue problem, rather than an ordinary 
differential equation based eigenvalue problem that is obtained in a classical analysis of primary instabilities 
in a shear flow that has strong gradients only along the wall-normal coordinate. The selection of grid and 
other aspects of the numerical solution were based on extensive experience with a similar class of flows (see, 
for instance. Refs. [64, 65]) and checks were made to ensure that the impact of variations with respect to those 
choices was negligible. Typically, the 2D eigenvalue problem corresponding to linear secondary instability 
was solved using 121 points along the wall-normal direction and between 33 to 65 grid points along the 
azimuthal direction. 



(a) Surface unrolling on a plane. The two dashed 
lines represent a meridian line on the cone and its 
imaginary print on the plane, respectively. 



(b) Correspondence between original and unrolled 
domains at a given X along with streamwise 
velocity contours. 


Fig. 2. Schematic of quasi-parallel secondary instability analysis. 


2. Results 

Computational results pertaining to both primary and secondary instabilities of the boundary layer flow over 
the cone are presented in this section. The evolution of stationary crossflow vortices excited by an 
azimuthally periodic array of localized roughness elements is described first in Section 2.1. The instability 
characteristics of selected vortices are investigated in Section 2.2. 

2 . 1 Evolution of Stationary Crossflow Vortices 

The stationary crossflow vortices leave an imprint on the cone surface in the form of longitudinal streaks in 
surface heat transfer as shown in Fig. 3(a). Immediately downstream of the forcing location, the vortex 
wavelength is approximately uniform along the azimuthal direction and corresponds to the forcing 
wavenumber of ^ = 60. The latter wavenumber was chosen to approximate the azimuthal wavelength of the 
crossflow modes found in the experiments of Craig and Saric [25]. Farther downstream from the forcing 


location, the streaks generated by the roughness element pattern progressively move towards the leeward ray 
and the streak spacing changes as a result of non-parallel flow effects. As discussed in Ref [66], there is a 
pronounced difference between the larger vortex wavelength near the sideline ray and the smaller wavelength 
on the leeward side. 

Axial velocity contours at selected streamwise stations are shown in Fig. 3(b). For the forcing location 
selected in this study (X= 0.18 m), the dominant crossflow vortices near the end of the computational domain 
lie between the sideline ray ((p = 90 deg) and the leeward ray ((p = 180 deg), analogous to that in the 
previously mentioned experimental investigations [21-25]. The velocity contours 3tX= 0.35 m and 0.40 m 
reveal the strongly nonlinear structure of crossflow vortices that qualitatively resembles the hot wire 
measurements of Craig and Saric [25]. 




(a) Contours of wall heat flux perturbation. (b) Axial velocity contours at selected locations. 

Fig. 3. Evolution of stationary crossflow vortices excited by azimuthally periodic array 
of localized roughness elements. 

Trajectories of selected vortices and the corresponding evolution of vortex amplitude along the cone length 
are shown in Figs. 4(a) and 4(b)-(c), respectively. Vortices 6 through 10 correspond to the highest velocity 
perturbations among those shown in Fig. 4(a). The amplitude plotted in Fig. 4(b) corresponds to the peak 
axial velocity perturbation over the cross section of the vortex, whereas the vortex amplitude shown in Fig. 
4(c) corresponds to the peak heat flux perturbation. Over most of the region shown in Fig. 4(b), vortex 8 has 
the highest velocity perturbation of all vortices. At X = 0.25 m, all three vortices have relatively small 
amplitudes, less than approximately two percent of the freestream velocity. By X= 0.30 m, the vortex 
amplitudes have increased to somewhere between flve to eight percent, whereas by X= 0.35 m, all three 
vortices have amplitudes in excess of 25 percent. The peak amplitude achieved by vortex 6 nearX= 0.35 m 
is almost 45 percent of the freestream velocity. The corresponding peak heat flux amplitudes in Fig. 4(c) are 
comparable to those measured in the experiments at Purdue University [21, 23, 67]. 
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(a) Selected vortex trajectories. 



(b) Vortex amplitudes measured by peak axial 
velocity perturbation. 



(c) Vortex amplitudes in terms of surface heat transfer. 

Fig. 4. Stationary crossflow vortices selected for secondary instability analysis. 

2. 2 Secondary Instability Characteristics 

The procedure for instability analysis outlined in Section 1 was applied to the basic state described in Section 
2.1. A representative spectrum of unstable modes is shown in Fig. 5, wherein the spatial growth rate of 
dominant instabilities has been plotted as a function of frequency for vortex 10 3tX= 0.25 m. AtX= 0.25 m, 
the amplitude of the primary crossflow vortex is relatively weak so that the predicted secondary instabilities 
are approximately equivalent to the primary instabilities of the unperturbed cone boundary layer with the 
additional constraint that the azimuthal wavenumber of the secondary mode be consistent with the local 
wavelength of the stationary vortex. 
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Fig. 5. Instability characteristics as a function of frequency at an upstream station (vortex 10, 

X=0.25 m). 

Figure 5 shows that the unstable secondary disturbances are concentrated in two separate frequency ranges. 
The low-frequency lobes (with peak frequencies below approximately 100 kHz) correspond to traveling 
crossflow instability modulated by the weak stationary crossflow vortex. Accordingly, they are labeled as 
TCI and TC2, respectively, where the numerical suffix identifies the index of the azimuthal harmonic 
corresponding to the stationary mode. Thus, the mode TCI corresponds to secondary modes with an 
azimuthal wavenumber equal to the local fundamental wavenumber of the stationary vortex, whereas mode 
TC2 corresponds to a secondary mode with an azimuthal wavenumber equal to twice the fundamental 
wavenumber. In a similar manner, the two high-frequency lobes correspond to second mode (i.e.. Mack 
mode) waves of the unperturbed boundary layer over the cone, and hence, have been labeled MMO and MMl , 
respectively. The mode MMO corresponds to axisymmetric Mack mode whereas mode MMl corresponds to 
an oblique Mack mode with an azimuthal wavenumber equal to that of the stationary crossflow vortex. 
Because the growth rate of second mode disturbances decreases rapidly with the wave angle, only the 
fundamental harmonic is unstable in this case. Representative mode shapes for the above families of 
disturbances are shown in Figs. 6(a) through 6(d). Peak velocity fluctuations associated with the Mack modes 
are relatively close to the surface, whereas those associated with the traveling crossflow modes are 
concentrated further away inside the boundary layer. In the remaining part of this paper, we focus mainly on 
the high-frequency disturbances since they appear to be linked more closely with the observed onset of 
laminar breakdown in the Purdue experiments [21, 22]. 

As the amplitude of the stationary crossflow vortex increases in the downstream direction, the modal topology 
gets progressively modified and it becomes difficult to attribute a clear mode type to each family of unstable 
secondary disturbances, especially to the low-frequency modes, unlike that in Fig. 5. Hence, the unstable 
modes at downstream locations are denoted simply using a numerical index, ordered approximately in the 
order of decreasing peak growth rate at the location of interest. A total of six dominant modes were selected 
for the plot in Fig. 7(a) 3tX= 0.30 m. Additional modes with relatively low peak growth rates have been 
omitted from the plot to avoid visual clutter. Representative mode shapes for the modes shown in Fig. 7(a) are 
plotted in Fig. 8. 
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Fig. 6. Representative mode shapes for secondary disturbances from Fig. 5. Coordinate Z 
denotes the azimuthal distance along the surface whereas Y denotes the wall-normal distance 
from the cone surface. 
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(a) 0.30 m. (b) 0.33 m. 

Fig. 7. Growth rate spectra for vortex 10 at selected downstream locations 

Results at intermediate locations (not shown herein) suggested that the maximum growth rate of the oblique 
Mack mode (mode MMl in Fig. 5) decreases with increasing X(and correspondingly, with a higher crossflow 
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vortex amplitude). Therefore, at locations sufficiently far downstream, the only dominant high-frequency 
mode related to the Mack mode instabilities at upstream stations corresponds to a modified version of mode 
MMO (corresponding to mode 1 in Fig. 7(a)). A comparison of Figs. 6(c) and 8(b) indicates how the original, 
azimuthally nearly uniform mode shape has been distorted by the stationary crossflow vortex. On the other 
hand, a new high frequency mode (mode 2 in Fig. 7(a)) that was not present at the upstream stations with low 
vortex amplitudes has emerged at 0.30 m. The peak growth rate of this new mode is small at this 
location, but becomes nearly the same as the growth rate of mode MMO farther downstream. The mode shape 
of this new mode bears a remarkable resemblance to the secondary instability mode driven by the spanwise 
shear of the crossflow vortex in swept-airfoil boundary layers [65, 68]. Thus, mode 2 apparently represents a 
genuine secondary instability of the finite amplitude crossflow vortex over the cone, and hence, is denoted 
alternatively by the nomenclature of SI. The peak-growth frequency of this mode is even higher than that of 
the mode MMO and the range of unstable frequencies 0.33 m is very broad, extending from almost 0 kHz 

to over 600 kHz. The large frequency bandwidth of this mode appears consistent with its interpretation as an 
instability of the relatively thin, inclined shear layer associated with the overturning contours of axial velocity 
associated with a strong crossflow vortex (Fig. 3(b)). 



(a) Mode TCI 



Z (mm) 

(b) Mode MMO 



(c) Mode SI (d) Mode 4 



Z (mm) 


(e) Mode 5 


Fig. 8. Representative |u’| mode shapes for secondary instabilities of vortex 10 at X= 0.30 m. 
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Since the unsteady disturbance spectra in Refs. [20-23] were measured via surface pressure transducers, the 
disturbance mode shapes for pressure fluctuation magnitudes for MMO and SI modes atX= 0.30 m are shown 
in Figs. 9(a) and 9(b), respectively. Analogous results forX= 0.33 m are plotted in Figs. 10(a) and 10(b). 
The flgures show that the peak pressure fluctuations for the MMO mode occur near the surface. As expected, 
azimuthal variations in surface pressure amplitudes becomes more pronounced as the vortex amplitude 
increases from 0.30 m toX= 0.33 m. The highest pressure fluctuations associated with the SI modes are 
concentrated away from the surface, approximately at the same location as the peak velocity fluctuations. 
The peak surface pressure fluctuations at the surface are approximately two times weaker than those at the 
interior peak. 



Z (mm) 



Z (itim) 


(a) Mode MMO. (b) Mode SI. 

Fig. 9. Ip’I mode shapes for secondary disturbance types MMO and SI (A^ = 0.3 m). 



Z (mm) Z (mm) 


(a) Mode MMO. (b) Mode SI. 

Fig. 10. Ip’I mode shapes for secondary disturbance types MMO and SI (X = 0.33 m). 

The cumulative growth characteristics for the most unstable MMO and SI modes are considered next. Figure 
11(a) shows the variation of the logarithmic amplification ratio, i.e., N-factor, along the cone length for 
selected frequencies of vortex 10 instabilities. Corresponding results for vortex 6 are shown in Fig. 1 1(b). 
Each red curve denotes a fixed frequency MMO mode, whereas the green curves denote the SI modes. For 
vortex 10, the most amplified MMO mode has a frequency of approximately 235 kHz, whereas the peak 
frequency of the SI mode is significantly higher, being equal to 320 kHz. The N-factor curves in Fig. 1 1(a) 
imply that the MMO modes of vortex 10 begin to amplify when the stationary crossflow modes are still very 
weak. Consistent with typical second (i.e.. Mack) mode behavior, the growth rates begin to decrease within a 
relatively short distance after their onset and then increase again as the amplitude of the stationary crossflow 
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vortex increases. For some disturbance frequencies, the latter increase is preceded by a modest decay in the 
amplitude of the MMO mode. In contrast to the MMO mode, the onset of SI modes occurs only after the 
amplitude of the crossflow vortex has become sufficiently large. In the range of approximately 0.35 m to 

X= 0.40 m (where vortex 10 has its largest amplitudes), the peak growth rates of the SI mode are larger than 
the peak growth rates of the MMO mode as indicated by the relative slopes of the corresponding N-factor 
curves. However, the peak N-factor of the SI modes at the end of the 0.40 meter length of the Purdue cone 
remains slightly smaller than the peak N-factor associated with the modified Mack modes. The N-factor 
differential between the two modes is approximately equal to one. The reader may also recall from Figs. 9 
and 10 that the mode MMO has a stronger signature at the surface than the SI mode. 

Compared to vortex 10, SI growth for vortex 6 begins farther downstream and the peak growth rates remain 
smaller. Consequently, the overall growth of the high-frequency disturbances over vortex 6 is entirely 
dominated by the MMO mode. The reasons behind the lower N-factors of SI modes of vortex 6 are not 
entirely clear. They might be related to the larger span wise wavelength of the crossflow vortices farther away 
from the leeward plane (since a larger wavelength would imply weaker azimuthal gradients of the basic state 
to facilitate energy transfer to the secondary instabilities) or the observation that the velocity amplitude of 
vortex 6 decreases rapidly after achieving its peak nearX= 0.37 m (Fig. 4(b)). The maximum Sl-mode N- 
factor for vortex 8 is in between that for vortices 6 and 10. 

The envelope of the N-factor curves for the MMO mode appears to level off between approximately 0.25 

andX= 0.3 m, and as the primary stationary vortex becomes even stronger downstream of this region, the N- 
factor envelope begins to rise again. In the absence of the stationary crossflow vortex, the maximum N-factors 
of second mode disturbances would remain below N = 2. The substantial increase in N-factors of the MMO 
modes in the region of high crossflow amplitudes illustrates the dramatic effect of the modulation of the MMO 
mode by the large amplitude crossflow vortex. 



X(m) 



0.4 

X(m) 


(a) Vortex 10. 


(b) Vortex 6. 


Fig. 11. N-factor evolution along selected vortices. 


Mode shape evolution for the MMO and SI modes of vortex 10 with the highest amplification ratios is shown 
in Fig. 12. The evolution of velocity disturbances associated with the MMO mode indicates how the peak 
fluctuations rise above the surface with increasing downstream distance. On the other hand, the mode shape 
of the SI mode remains very similar throughout the evolution of this mode, being concentrated along the 
inclined shear layer as mentioned previously. 
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(a) Mode MMO,X= 0.25 m 




(b) Mode MMO, 0.28 m (c) Mode SI, X= 0.30 m 
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(d) Mode MMO, X= 0.325 m 


(e) Mode SI, X= 0.325 m 


Fig. 12. Streamwise evolution of |u’| mode shapes for dominant instabilities of vortex 10 

(continued). 
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(f) Mode MMO, X= 0.35 m (g) Mode SI, X= 0.35 m 



> 



2 (mm) 

(h) Mode MMO, 0.375 m 


(i) Mode SI, 0.375 m 


Fig. 12. Concluded. 


A comparison of the results plotted in Figs. 1 1(a) and 1 1(b) shows that the peak frequency of MMO mode 
decreases from vortex 6 to vortex 10 because of the increase in boundary layer thickness from windward to 
leeward planes. On the other hand, the frequency of the SI Mode is strongly influenced by the characteristics 
of the inclined shear layer, which are determined by the wavelength and amplitude of the stationary crossflow 
vortex. Thus, one finds, in fact, that the peak disturbance frequency of the SI mode as predicted via the quasi- 
parallel stability theory employed herein increases from vortex 6 to vortex 10. 

Both experiments [69] and computations [70-75] of crossflow breakdown in subsonic swept-wing boundary 
layers have shown that the high-frequency secondary instabilities can rapidly amplify by orders of magnitude 
to initiate the onset of transition. Even though the SI modes have large growth rates in the present, hypersonic 
case, the modified Mack modes also play an important role and may in fact correspond to higher 
amplification ratios than the SI modes. Thus, future computations of nonlinear interactions involving the 
broader set of high-frequency modes would be of great interest. 

Figure 13 shows the azimuthal variation in growth rate spectra for modes MMO atX= 0.33 m. In particular, 
growth rate spectra for vortices 6, 8, and 10 are shown as functions of the disturbance frequency. The peak 
frequency of the MMO modes is the highest for vortex 6 and lowest for vortex 10. The variation in peak 
frequency is qualitatively consistent with the accompanying increase in underlying boundary layer thickness. 
Unlike the MMO modes, the peak frequency of the SI modes does not show any pronounced variation from 
vortex 6 to vortex 8. On the other hand, unlike the MMO modes, the peak growth rate of the SI Modes 
increases substantially from vortex 6 to vortex 10, in spite of the fact that the peak velocity amplitude of 
vortex 10 at X=0. 33m is lower than the amplitudes of the other two vortices. The variation in growth rate 
may be related to the smaller azimuthal wavelength of vortex 10, which implies potentially stronger azimuthal 
gradients in the basic state that drive this particular mode of secondary instability. 
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Fig. 13. Variation in growth rate spectra along azimuthal direction (X = 0.33 m). The curves 
denoted as V06, VOS, and V10 indicate the growth rate variations for vortices 6, 8, and 10, 
respectively. 


3. Summary and Concluding Remarks 

Quiet tunnel experiments at Purdue University [21] have detected high-frequency disturbances in the vicinity 
of an apparent breakdown of large amplitude stationary crossflow vortices during their experiments on a 
yawed circular cone in the Mach 6 B AMQT facility at Purdue University. The observed dependence of these 
high-frequency disturbances on a threshold stationary wave amplitude and their localized azimuthal signature 
led these investigators to suggest that these disturbances were caused by the secondary instability of the 
stationary crossflow modes. In contrast, Munoz et al. [20] attributed the high-frequency measurements in a 
conventional facility to second mode waves. The instability analysis in this paper shows that the high- 
frequency disturbances may be attributed to both second mode waves modified by finite amplitude crossflow 
vortices and a genuine secondary instability of the stationary vortices. 

At sufficiently high stationary crossflow amplitudes, the original second mode (i.e.. Mack mode) disturbances 
are further destabilized, resulting in peak growth rates that are significantly larger than those of second modes 
in the absence of any crossflow vortices. Large amplitude crossflow vortices also lead to a genuinely new 
instability mode (mode SI) that corresponds to secondary instability of the azimuthally modulated boundary 
layer. Unlike the second mode disturbances, the SI mode is unstable over a broad range of frequencies and its 
mode shapes are similar to the secondary instabilities driven by spanwise basic-state shear in subsonic swept 
wing boundary layers. The frequencies of the most amplified modified Mack modes and the SI mode are in 
the range of 200 kHz to 300 kHz, with noticeable variation depending on the location of interest and the 
amplitude of the crossflow vortex. The above frequency range is approximately consistent with the 
frequencies measured in the experiments at Purdue University [21-23] and the Technical University of 
Braunschweig [20]. However, additional computations to assess the role of the modeling approximations 
employed in the present work are deemed highly desirable. In particular, the effect of non-parallel evolution 
of the basic state for secondary instability must be examined. Given the richness of the unstable secondary 
modes within a relatively small range of phase speeds, it is possible that the non-parallel effects would lead to 
significant changes in the modal topology and, hence, in the evolution of the N-factor envelopes. 
Computations of this type are currently in progress and the results thereof will be reported in a separate paper. 

Finally, the spatial structure of the computed SI mode is similar to the measured mode shape of the so-called 
high-frequency modes in the quite tunnel experiments [25] at Texas A&M University. However, the present 
SI mode has substantially higher frequencies than the 80-120 kHz range measured in those experiments. 
Thus, the origin of the latter modes remains to be determined. 
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The presence of multiple SI modes in the Mach 6 configuration of interest is analogous to the secondary 
instability of crossflow vortices over subsonic and supersonic airfoils. Computations for the subsonic 
boundary layers have shown that nonlinear effects tend to reduce the growth of the secondary instability 
modes. However, if the amplitudes become sufficiently large, then a rapid breakdown to turbulence can 
ensue. Whether or not a similar scenario prevails at hypersonic Mach numbers, where the modified Mack 
modes also play an important role, remains to be seen. The follow-on work will focus on the breakdown 
mechanisms and the strongly nonlinear phase of transition to turbulence at hypersonic speeds. 
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